Heteroclinic travelling waves in convex FPU-type chains 

Michael Herrmann* Jens D.M. Rademacher^ 

November 3, 2009 

Abstract 

We consider infinite FPU-type atomic chains with general convex potentials and study the exis- 
tence of monotone fronts that are heteroclinic travelling waves connecting constant asymptotic states, 
looss showed that small amplitude fronts bifurcate from convex-concave turning points of the force. 
In this paper we prove that fronts exist for any asymptotic states that satisfy certain constraints. 
For potentials whose derivative has exactly one turning point these constraints precisely mean that 
the front corresponds to an energy conserving supersonic shock of the 'p-system', which is the naive 
hyperbolic continuum limit of the chain. The proof goes via minimizing an action functional for 
the deviation from this discontinuous shock profile. We also discuss qualitative properties and the 
numerical computation of fronts. 
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1 Introduction 

We consider infinite chains of identical particles as plotted in Figure [T] These are nearest neighbour 
coupled in a convex potential $ : R — > R by Newton's equations 

Xa = ^'{Xa+l - Xa) - ^'{Xa - Xa-l), (1) 

where ' = ^ is the time derivative, Xa{t) the atomic position, a G Z the particle index. 

AAM — • — ^/\A/^ — • — ^/\A^ 



Figure 1: The atomic chain with nearest neighbour interaction. 



Such chains model particles connected by springs in one dimension and serve as simplified models for 
crystals and solids. In their seminal paper [FPU55| Fermi, Pasta and Ulam studied such chains assuming 
that the interaction potential 4> contains at most quartic terms. We consider convex $ with nonlinear 
force function $' and allow for turning points, i.e., points where <!>'" = 0, but still refer to ([T]) as FPU 
chains. In fact, the existence of fronts, which is studied in this paper, requires that has at least one 
turning point, see |IIR08| . and this excludes, for instance, the famous Toda potential. 

One focus of the research concerns solutions with coherent structures, in particular travelling waves, 
for which there exists a smooth profile that travels with constant speed and shape through the chain. 
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The main types of travelling waves are periodic wave trains, homoclinic solitons (or solitary waves), and 
heteroclinic fronts. Mathematically, rigorous existence proofs of such solutions are an important basic 
issue. During the last two decades a lot of research addressed the existence of solitons and wave trains: 
[FW94[ IFV991 IHer0 5 , DHM06| establish the existence of such waves by solving constrained optimisation 
problems, |SW97[ iPPOOi iPanOSi ISZ07j apply the Mountain Path Theorem to the action integral for 
travelling waves, and |IooOO| uses center manifold reduction with respect to the spatial dynamics. 

In comparison, little is known rigorously for fronts. For (non-smooth and non-convex) double-well 
potentials composed of the same quadratic parabolas, the existence of fronts connecting oscillatory states 
has been recently shown in |SZ08j . Such fronts can be interpreted as phase transitions and more physical 
results can be found for instance in |SCC05[ ITVOSj . 

In these cases the connection between fronts and shocks in the naive continuum limit of ([T]) was 
crucial. This limit is the so-called p-system formed by the hyperbolic conservation laws for mass and 
momentum. Shocks in the p-system come in different types given by the relation of their speed and 
the sound speed of the asymptotic states: shocks that are faster (slower) than these sound speeds are 
called supersonic (subsonic). The fronts found in [SZ08| for the quadratic double- well case correspond to 
subsonic shocks when taking the average of the asymptotic oscillations. It turns out that these shocks 
also balance the energy jump relation, and we refer to such shocks as conservative. 

The connection to p-system shocks is also crucial for our result, but we solely consider convex potentials 
and the fronts we find are supersonic, conservative, and monotone with constant asymptotic states, 
see Figure [2J For such potentials, the only previous result concerning fronts we are aware of is the 
bifurcation result by looss in |IooOO| for supersonic fronts of small amplitude connecting constant states 
near a convex-concave turning point of A physical interpretation of these solutions is a deformation 
of the solid or crystal, though the fronts can travel in either direction for the same asymptotic states. 
Related 'pulsating' travelling waves have been studied in [IJ05j . Here the profile of the wave changes 
periodically while moving and the asymptotic states are wave trains. It was shown that in the truncated 
normal form on the centre manifold, pulsating front solutions bifurcate from convex-concave turning 
points of 

The analytical investigations in this paper are motivated by the numerical simulations of atomistic 
Riemann problems which have recently been studied by the authors in |HR08| . We observed fronts in 
numerical simulations of ([1]) even for initial data that are far from the data of a front, see Figure[21 Hence, 
fronts are dynamically stable and provide fundamental building blocks for atomistic Riemann solvers. 
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Figure 2: Fronts appear in the numerical simulations of FPU chains: snapshots of the atomic data with 
monotone front around particle a = 1613. 



The main result of this paper is the following. 

Theorem 1.1. For all convex and smooth potentials $ the following assertions are satisfied: 

1. Each front in the chain corresponds to a conservative shock in the p-system. 

2. For each supersonic conservative shock in the p-system that satisfies a certain area condition there 
exists a corresponding monotone front in the chain. 



Here 'correspond' means that shock and front have the same wave speed and asymptotic states. 
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The precise assertions of this theorem depend on the set of turning points of <i>': 

1. If $' has no turning points, then all shocks arc non-conscrvativc and fronts do not exist. 

2. For strictly convex-concave <&' all conservative shocks are supersonic and satisfy the area condition. 
In this case Theorem 11.11 implies that each conservative shock can be realised by an atomistic front. 
In particular, this allows for arbitrary large jumps between the asymptotic states and extends the 
bifurcation result from |IooOO| . 

3. For strictly concave-convex <&' all conservative shocks are subsonic and Theorem 1 1 . 1 1 does not imply 
an existence result for fronts. 

4. If $' has more than one turning point, then Theorem 11.11 provides the existence of fronts for a 
proper subset of all conservative shocks. 

We emphasise that Theorem 11.11 does not imply the uniqueness of fronts with prescribed asymptotic 
states. Moreover, not every conservative shocks in the p-system corresponds to a front in the FPU 
chain. For functions with one turning point the simulations in |HR08| rather suggest that only 
supersonic shocks allow for an atomistic realization as a (dynamically stable) front. This is supported 
by the bifurcation result in [IooOO| , which in fact disproves the existence of subsonic conservative shocks 
with small jump heights, and our results, which imply that subsonic shocks cannot minimize the action, 
compare the discussion at the end of 321 

We also mention the results about atomistic fronts by Aubry and Proville. In |AP07| they study FPU 
chains with concave $' - so all shocks in the p-systcm are non-conservative - and allow for atomistic 
fronts by adding some damping force that respects uniform motion. By numerical simulations they found 
oscillations in the tails of such fronts and observed that these oscillations persist when the damping tends 
to zero. We therefore expect that non-conservative shock data can be related to heteroclinic connections 
of wave trains. 

The proof of the first part of Theorem 11.11 is fairly straightforward and contained in Lemma 12.11 below. 
Proving the second part, which is mainly contained in Theorem 13.141 requires a number of prepara- 
tory steps, and uses the convexity of (f>, the supersonic front speed, and the area condition in various 
fundamental steps. To guide the reader we give an overview of the key ideas and steps for the proof. 

1. We use the a priori knowledge of the front speed in order to reformulate the problem as a nonlinear 
fixed point equation for a suitably normalised profile. 

2. We identify an action functional for the deviation from the discontinuous shock profile such that 
the fixed point equation is the corresponding Euler-Lagrange equation. 

3. We use the invariancc of the set of monotone profiles under the gradient flow of the Lagrangian to 
connect stationary points in this set with fronts. 

4. We establish bounds for the action functional and use the direct approach to show that the La- 
grangian attains its minimum in the set of monotone profiles. 

This paper is organised as follows. In ^ we give a more formal background for fronts and shocks of 
the p-system, and contains the proof of the existence result. In iJHwc prove that these fronts decay 
exponentially to their asymptotic states; the proof also applies to supersonic solitons with monotone 
tails and extends the known results there. In addition, we discuss numerical computations based on 
the gradient flow, and hereby illustrate the role of area condition, supersonic front speed, and energy 
conservation. 
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2 Fronts, shocks and normalisation 



For our purposes it is convenient to use the atomic distances Tq, — Xa+i — Xa and velocities Va = ia as 
the basic variables, changing ([T|) to the system 

ra = Va+l - Va , Va = ^'(''q) - $'(ra-l). 

Travelling waves are exact solutions to the infinite chain ([T|) that depend on a single phase variable 
if = a — at, where the speed a equals the phase velocity. In terms of atomic distances and velocities 
travelling waves can be written as 

rait)^R{ip), Vait)^V{ip), 

where the profile functions R and V solve the advance-delay differential equations 

a^R{^) + V{^ + l)-V{^)^Q, + -<i>'(i?(^-l)) =0, (2) 

which imply the energy law 

^^{kV^iv) + miv))) + ^'{R{^))Vi^ + 1) - <I>'(i?(^ - l))V{^) = 0. (3) 

The travelling waves in our context are fronts, which connect two different constant asymptotic states 
for if —>■ —oo and ip oo, i.e., 

R,VeOnL°° with R{ip) r± and V{ip) v±, (4) 

where and C'^ denote the usual function spaces on the real line. The other variants of travelling waves 
mentioned in the introduction are wave trains, with periodic R and V, and solitons, that are localised 
over a constant background state, i.e., r_|_ = r_ and — w_. 

In |IooOO| it has been proven that fronts bifurcate from turning points of with <^^^\r^) ^ 
if and only if $'^'(r») < 0. These fronts have strictly monotone profiles that converge exponentially to 
their asymptotic states. Our main existence result extends this in the sense that we prove existence of 
fronts for any left and right states that satisfy certain constraints and can have large amplitude. 

As a first result we establish necessary conditions for the existence of fronts, which reflect that each front 
converges to a shock in a large-scale limit. These conditions turn out to be the jump conditions for an 
energy conserving shock of the p-system (see below) , and show that the speed cr of a front is completely 
determined by its asymptotic states. 

For any atomic observable ij: = '0(r, v) we define the mean and the jump by 

{i^) ■=\{'^{r+,v+) + ^{r_,V-)) and m := ^:{r+, v+) - ^(r_, v.), 

respectively, where r± and v± denote the asymptotic values of a front, see (^1). 

Lemma 2.1. Suppose ttiat R and V as in ^ solve the front equation ([2]). Then, 

a|[r]| + IHI = 0, a|H| + |[<f>'(r)]=0, + + 0. (5) 

Proof. Integrate over a finite interval [—n^ n] gives 

n-t-l n — n+1 n+1 

<j{R{n) - R{-n)) = - ( V {ip) Aip + [ V [ip) d^p = [ V{^)d^~ [ V^M d(^, 



-n+1 



and the limit n — > oo yields ([5])i. The jump conditions ^2 and (Oa follow similarly from ([2])2 and 
©. □ 

We mention that the jump conditions ^ are also derived in |AP07j and that a similar method is used 
in |Ser07| to derive jump conditions from purely discrete front equations. 
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2.1 The p-system and conservative shocks 

The p-system can be formally derived from the FPU chain as a continuum limit in the hyperbolic scaling 
of the microscopic coordinates t and a. This scaling bridges to the macroscopic time t = et and particle 
index a = ea, where < e <C 1 is the scaling parameter. It is natural to scale the atomic positions in 
the same way, i.e. x = ex, which leaves atomic distances and velocities scale invariant. 

Substituting an ansatz of macroscopic fields r(t,a) and v(t,a) such that ra{t) = r{et, ea), Va(t) = 
v[et, ea) into ([1]) and taking the limit e — > yields the macroscopic conservation laws for mass and 
momentum 

djr^d-sv^Q, djv - d^^' {r) = 0, (6) 

which is the aforementioned p-system. It is well known that the p-system is hyperbolic if and only if <& 
is convex, and that for smooth solutions the energy is conserved via 

dj + $(r)) - {v $'(r)) = 0. 

Next we summarise some facts about shocks in the p-system. For more details we refer to standard 
textbooks |Smo94i [DafOOl ILeF02] . and to |HR08| for conservative shocks. 

A shock connecting a left state u- — (r_ , f _ ) to a right state uj^ ~ {r+ , ij+ ) propagates with a 
constant shock speed a so that u- and m+ satisfy the Rankine-Hugeniot jump conditions for mass and 
momentum ([5|)i and Corresponding to the two characteristic velocities of ([6]) one distinguishes 

between 1-shocks with cr < and 2-shocks with cr > 0. 

An important observation is that for either convex or concave flux the jump conditions and 
^2 imply that the jump condition for the energy must be violated, that means 

<y\[\v^ + Hr)]\ + b^\r)]\^Q. (7) 

For non-convex or non-concave <&', however, it is possible to find conservative shocks that satisfy all three 
jump conditions from (O, see Theorem 5.1 in |HR08j . Eliminating \v\ and a in ([5]) shows that each 
conservative shock satisfies 

J(r_,r+):=|[$(r)]-|H(<I>'(r))=0. (8) 

Conversely, for each solution to ([8]) there exist, up to Galilean transformations, exactly two corresponding 
conservative shocks which differ only in sgn|u] = — sgncr and propagate in opposite directions. The 
geometric interpretation of jS]) is that the signed area between the graphs of and the secant line 
through $'(r_) and <I>'(r+) is zero over [r_,r+], compare Figure [4] below. 

For the harmonic potential (f>(r) = the p-system is linearly degenerate, hence all shocks are 
conservative. In the nonlinear case, however, the solution set to is either empty or the union of curves 
in the (r_, r+)-plane, see Figure[3]for an example. In particular, precisely from turning points of there 
bifurcate curves of conservative shocks, see |HR08j . 

Shocks come in different types determined by the relation of the shock speed a and the sound speeds 
A(r) := ±yj^"{r) evaluated in r_ and r+, where the signs of both the sound speed and the shock speed 
are taken negative (positive) when considering 1-shocks (2-shocks). More precisely, a shock connecting 
u- to u+ with speed a is called 

1. compressive, or Lax shock, if A(r_) > cr > A(r+), 

2. rarefaction shock, if A(r_) < a < \{r+), 

3. supersonic, or fast undercompressive, if \a\ > |A(r_)| and \a\ > |A(r+)|, 

4. subsonic, or slow undercompressive, if |(t| < |A(r_)| and \a\ < \X{r^)\, 

5. sonic, if \cr\ = |A(r_)| or \a\ = |A(r+)|. 

For all sufficiently small jump heights |r+ — r_ | + \v+ — \ one can show that compressive and rarefaction 
shocks have negative and positive, respectively, energy production in ([7]). Consequently, all conservative 
shocks with small jump height are either sonic or undercompressive. The conservative shocks bifurcating 
from turning points of $' are supersonic if $'-^-'(r*) < and subsonic if $^^^(r*) > 0, see jHR08] . 
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Figure 3: Conservative shocks for the potential <i>(r + 1) = r + ^r^ + — \ cos(2r) + sin(3r): sohd 
hnes form the sohition set to J^(r_,r+) = 0, i.e., the distance data for conservative shocks; bullets lie at 
local extrema of the curves in the coordinate directions, and are the points where the shock type changes. 



2.2 Normalisation and reformulation 



Next we show that, after a suitable renormalisation, the front equations for R and V are equivalent to 
a nonlinear fixed point equation for a normalised profile W with asymptotic states ±1. To this end we 
introduce the averaging operator 



urn d^, 



which satisfies 



dip 



u( 



1 

V-2 



lim {AU^iip) = lim U{ip). 



C^— >±oo 



ip — >-±oo 



(9) 



For given asymptotic front states (r±, v±) we further introduce a normalised potential $ by 



W{r) 



^{r) + \\[r]\w 



|[$'(r)]| 



which is convex on [—1, 1] and strictly convex if $ is. Moreover, the jump conditions from ([5]) imply 
$'(±1) = ±1 and $(1) = $(—1), and the shock corresponding to (r±, v±) is supersonic if and only if 
$"(1) < 1 and $"(-1) < 1. 

Lemma 2.2. For fixed asymptotic states (r±, v±) there is a one-to-one correspondence between front 
solutions to ^ and solutions W to 



W = A<^>'{AW), 

with W{ip) ±1. 

Proof. Suppose that R and V solve ([2]). Then there exists a normalised profile W such that 

V ^{v) + \\[v]\W. 
By ([9]) the first equation from ([2]) is equivalent to 

-cjR = AV + c^{v) + \ [v]\AW + c 
with i?(iy9) = R{(p — ^) and an integration constant c. Passing to ^ — oo and — > oo we find 



(10) 



(11) 
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and this gives 

R= {r) + ^\[r]\AW. (12) 
The second equation from ([2|) transforms into 

-aV ^ A<i>'{R) +d, (13) 
and inspecting the asymptotic values we find 

-a\[v]\^\[<^'{r)l d=-a{v)-{^'ir)). (14) 
Inserting ^ and ^ into ^ gives 

W = + ^\[r]\AW) - ($'(r))), 

which is (fTU|) due to P^ . FinaUy, it is easy to see that each solution to (fTU)) determines a front via pT|) 
and ([H]). □ 

The equivalence between the lattice equation Q for travelling waves and the nonlinear eigenvalue 
problem a-^W = A^'{AW) is well established for wave trains and solitons, compare for instance |FW94[ 
IFV99| ■ but in the case of fronts we have not seen this formulation before. An important difference is that 
for fronts the speed a is completely determined by the asymptotic states and can hence be removed from 
the problem. 

We conclude with some remarks concerning the regularity of fronts. Suppose that W G L°° solves ^TU\\ . 
Using ^ , compare also Lemma 13.31 below, and the smoothness of $ we find 



d<p'- 



■W = \7[d>"{AW)VW], (15) 



and hence W £ W^'°° C C. Moreover, we infer that W & and, differentiating further, that W is 
as smooth as $. Finally, exploiting pT|) and ([12]) we arrive at the following result and refer to for 
further qualitative properties of fronts. 

Remark 2.3. Each solution W G L°° to (jlOp is contained in and provides a smooth solution to (|2|) 
with at least R,V £ O r\L°°. 

3 Proof of the existence result 

In order to prove the existence of monotone solutions W to the fixed point problem (fTU|) we consider only 
normalised potentials <E> = $ and rely on the following standing assumption. This involves the function 

g^{w)-~ J V- $'(w) dv = $(-1) - $(u;) + ^w^ - i, (16) 
-1 

which measures the signed area between the identity and the graph of see Figured 
Assumption 3.1. The potential $ satisfies 

(R) regularity: $ is twice continuously differentiable on [—1, 1], 

(N) normalisation: $'(-1) = -1 and $'(1) = 1, 

(C) convexity: $" is nonnegative on [—1, 1], 
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(E) conservation of energy: $(—1) = $(1), 

(S) supersonic front speed: $"( — 1) < 1 and $"(1) < 1. 

(A) positive signed area: g^{w) > for all w G (—1, !)• 

We mention that (N) and (R) are made for convenience, (E) is necessary according to Lemma 12.11 
and (C) is needed for the monotonicity of fronts. The status of (S) and (A) is more dehcate: If $' has 
only one turning point, then (S) coincides with the bifurcation condition in [fooOO] and imphes (A), see 
Remark 13.21 below. If $' has more than one turning point, then (A) is independent of (S). In any case 
both (S) and (A) are closely related to the existence of action-minimising fronts as discussed in more 
detail at the end of fJS] and illustrated in i|4.2l 




Figure 4: Two examples for force functions $' that satisfy Assumption [331 the first one is prototypical for 
the supersonic fronts with small amplitudes that bifurcate from convex-concave turning points. Condition 
(A) precisely means that for each — 1 < < 1 the signed area between the identity and the graph of $' is 
positive in the stripe [—1, w]; the signed area vanishes in the stripe [—1, 1] if and only if the corresponding 
p-system shock conserves the energy via (E). 



Recall that turning points of $' are exactly those points in which <!>" changes its monotonicity. More 
precisely, convex-concave (concave-convex) turning points correspond to local maximisers (minimiser) of 
We now summarize the main observations concerning the turning points of 

Remark 3.2. Suppose that (N), (R), and (S) are satisfied. Then $' has at least one convex-concave 
turning point in (—1, 1). Moreover, if ^' has no further turning point in (—1, 1) then (E) implies (A). 

Proof. Suppose for contradiction that max <I>"|[_x_i] < 1. Then, for all w G [—1, 1] we find g%{w) = 
1 - > and hence = 5$(-l) < .g<i.'(w) < .9*'(+l) = 0. This implies g'^{w) = for all 

w G [—1, 1], and hence a contradiction to (S). Therefore, possesses a local maximiser w G (—1, 1), 
which is a convex-concave turning point of 

Towards (A) suppose additionally (E) and that is convex-concave in [—1, 1]. Then, 5$(1) = 
(?$(— 1) = 0, and the monotonicity properties of imply that there exist —l<wi<W2<i such that 
^"{w) < 1 and > 1 for all w £ Ii ~ [—1, wi] U [w2, 1] and w e I2 = [wi, W2], respectively. In 

particular, g^ is convex in Ii but concave in I2, and combining this with = ^$(±1) = .g^(±l) and 
g'^i±l) > 0, we readily find g<s>iw) > for all w G (-1, +1). □ 

3.1 The action functional 

In this section we cast the problem in a variational framework and show that 

W = T[W], T[W] ~ A^'{AW), 
is the Euler-Lagrange equation for an action integral. To this end we introduce the afhne Banach space 
n := (functions : M ^ K such that = Wsh - e lH, 
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where the reference profile Wsh{^) ~ sgn{(p) is the shock profile that connects —1 to 1. We also define 
the Lagrangian on Ti. by 

C{W) = J L{W, AW) - L(W^,h, ^M4h) dip, L{v, r) = ^v^ - $(r), (17) 

and a formal calculation shows = W — T[W]. Notice that the ansatz W € L'^ arises naturally 

when computing the asymptotic behaviour for ip —^ ±oo. In fact, since AW decays like W and due to 
$'(±1) = -1 the integrand in ^ behaves like i(l - {±1))W (ipf . 

In the remainder of this section we show that both £ and T are well-defined on H. We start with some 
properties of the averaging operator A that are proven in |Her09j . 

Lemma 3.3. The linear operator A is well-defined on L^, 1 < p < oo, and has the following properties: 

1. A maps into VJ^'P with (AWYip) W{p+ 1/2) - W{p- 1/2), 

2. A maps \? into \? n L°° with \\A^\\oo<\W\\2 and |i2<|i II2, 

3. A is self-adjoint on L^, 

4-. If a sequence {Wn)n converges weakly in to some limit Woo, then {AWn)n converges strongly 
and pointwise to AWoo in L^([— M, M]) for each M < 00. 

Our next result shows T[W] - W E L'^ for all W e H and implies that T : H ^ H is well-defined. 

Lemma 3.4. For any W £ H we have 

W-^'iAW) e L^, W~A'1>'{AW) e L^, (18) 

as well as 

J W^l~^'iAW)A^Idip^ J W^l- A<P'iAW)fidip (19) 

for all test functions fj, G . 
Proof With 

\<S>'{{AW)i^)) ^^'{(AWsi,){ip))\<C\AW{ip)l C = sup{|$"(u.)| : \w\ < max{l, \\ AW \\oo}} 
for all ip, and ^'(AWshif)) = Wshif) for \ip\ > \, we find (with the usual cut-off function x) 
\W --<^'{AW)\ <W + x,i i\W^h-'^'{AWsY,)\ + C\AW\. 



2 ' 2J 



This gives (|18P i since W , AW G L^, and the proof of is similar. Now let M > be arbitrary. Since 
A is symmetric, see Lemma 13.31 we write 



M M 



J <i>'iAW)Afidip = J fiA{x[-M.M]^'{AW))dip^ J fiA1>'{AW)dip + EM, (20) 



-M R -M 

where the error term satisfies 

\Em\<C J \ti\d.p+ J \fi\d^ 



C = sup {!$'(«;) I : \w\ <\\AW\\oo }• 



Finally, we add J^^j Wfi dip on both sides of ([20]) , and taking the limit ilf ^ 00 yields (fT9|) . □ 
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We now prove that L is well-defined and Gateaux-difFerentiable on Ji. 
Remark 3.5. For each Wi^ the pointwise normalised potential 

$^oH := ^{wq + w)- $'(wo)w - <^{wo) 
satisfies < ^wg{w) < ^liP' with c = sup$"|[j^,g 

Lemma 3.6. The functional C is well-defined on Ti. Moreover, it is Gateaux differentiable with locally 
Lipschitz continuous derivative dC\W] = W — T[W]. 

Proof. By construction, we have 

L{W, AW) - L{Wsb, AWsh) = - W,^W + <i>'{AW,^)AW - -^awA^W) 

with pointwise normalised potentials ^aw^^, as in Remark [531 Using ([19)) with W — Wsh and ii = W we 
find 



/ 



W,hW ~ ^\AW,i,)AW d^ = j (w^sh - A^'{AW,iS) W dip, 



with well-defined integrals on both sides. Using C :~ sup$"|[_i ij < oo wc estimate 
J ^W^dip < oo, J <pAwAAW)dip < C \\AW\\l< C \\W\\l, 

R R 

and conclude that C{W) is wcU-dcfincd for all W £ H. Moreover, (fT^ implies W — T[W] & L^, and in 
view of (|19p we readily verify the formula for dC as well as the claimed Lipschitz property. □ 

3.2 Monotone and heteroclinic functions 

In order to introduce a refined ansatz space for the front profile W we start with some preliminary 
remarks concerning the weak formulation of monotonicity, asymptotic values, and L°° bounds. 

Let [ipo, oo) be some interval, U G L°°([(^07 oo)) be fixed, and let ^ > denote an arbitrary smooth 
test function with compact support in [ipo, oo). Then, the function U is said to 

1. be increasing, if {U{(p + (p) — U{ip))^{(p) dip > for all shifts (p > and all fj,, 

2. be nonnegative, if J^U {(p) fi{ip) dip > for all /i, 

3. take values in [ui, U2], if U2 > U{ip)fi{ip) dip > ui for all > with ^{ip) dp = 1, 

4. have the asymptotic value it-|_ for p — > 00, if J^U{ip + ip)^{ip) dip Jjg^fi{p)dp as <^ — > 00 for 
all ^. 

Note that all these properties are preserved under weak convergence, and that similar statements hold 
on (00, po]- 

We aim at establishing the existence of fronts in the convex set 

C := e ?Y is increasing with W{p) ±1 as ip ^ ±00, | 
where the convergence conditions mean that W has asymptotic values —1 and +1 in the above sense. 

A particular problem in our subsequent analysis is to control the relative shift between a profile function 
W and the reference profile Wsh- For this reason wc introduce the set of all pinned profile functions by 



Co :={weC with ;^!["'+-) "^1^^^ 

I I (-00,0] takes values m [—1, OJ 
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Obviously, for each W € C there exists a phase shift (po such that W{- — (po) G Co, and this phase shift is 
unique if W is strictly increasing. Moreover, we have G Co if and only if e Co with 

Co •= e such that ^^1°'+°°) decreasing with values in [0, +1] 1 
] W^l(-oo,o] is decreasing with values in [—1, 0] J 



Notice that the set Co is convex and closed under weak convergence in L 
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3.3 Variational approach 

The starting point for our variational existence result is the invariance of C under the action of T. In fact, 
it is straightforward that the averaging operator A respects both the monotonicity and asymptotic values 
of W. In addition, since $' is a monotone map on [—1, 1] with $'(±1) = ±1 the set C is also invariant 
under the pointwise application of Notice, however, that Co is not invariant under the action of T 
since for (continuous) W with W{0) = in general we have T[VK](0) ^ 0. 

Next we consider the gradient flow of C inH, that is 

-d.w = -dC[W] = ~W + T[W] (21) 

with flow time s. According to Lemma[321 the initial value problem for this 7-^-valued ODE is well-defined, 
and we readily verify that 

■£C{W) = {dC{W), -^W) =-\\W~ T[W] \\l . (22) 
In particular, each stationary point of pip is a fixed point of T. 

The key observation is that C is an invariant set for (j21|) . To sec this we introduce the corresponding 
Euler scheme with step size < A < 1, that is 

W ^ W - XdC[W] = (1 - X)W + XT[W]. (23) 

This scheme leaves C invariant because W ^ C implies T[W] S C, and since C is convex and closed in TC 
the claim follows with A ^ 0. As a consequence we obtain the following result. 

Lemma 3.7. Each local minimiser of C in C is a solution to the front equation W = T[W]. 

We proceed with two remarks: (i) The front equation can be considered as the Euler-Lagrange 
equation for C corresponding to variations W W + SW with arbitrary W € i.^. Therefore, the fact 
that minimiscrs of £ in C solve the front equation is not clear a priori but follows from of the invariance 
properties of C and the dissipation inequality (|22|) . (m) There is no analogue to Lemma [3.71 for local 
maximiscrs. This follows from the next Remark and reflects that C is not invariant under the negative 
gradient flow of C 

Remark 3.8. The shock profile M4h is a local maximiser for C inC but does not satisfy the front equation. 

Proof. Let G C be arbitrary with W ^ Wsh- For < e < 1 we define We = (1 — £)Wsh -\- eW and note 
that We G C. Lemma 1321 yields 

A —C{We)\e=o = (a£[I^.h], W-Wsh)^~l {Wsh - T[Wsh]){Wsh - W) dp. 

M 

On the one hand, we have T[Wsh] G C with W^hif) < TlW^hlif) for -1 < </? < 0, Wsh{^) > T[Wsh](</') 
for < < 1, and M4h(¥') = T[M4h](<^) for \ip\ > 1. On the other hand, 1^ G C implies Wshi^p) < W{ip) 
and Wship) > W{(p) for (p < and f > 0, respectively. From this we conclude A < and W ^ Wsh gives 
A < 0. Finally, Wsh does not equal T[Wsh] as the latter is continuous. □ 

Our strategy for proving the existence of fronts is to show that C attains its minimum in C. To this 
end we follow the direct approach, and prove that minimising sequences are precompact. 
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3.4 Bounds for L 

Towards compactness results for minimising sequences for £ in C, we first prove that corresponding 
sequences in Co also minimise £. To this end we start with a technical result. 

Lemma 3.9. Let U £ L°° be monotone with asymptotic values u_ and for if — > —oo and (f — > oo, 
respectively. Then we have 

U{ip + ip) - U{(p)difi = {u+ - u-)'p. (24) 



R 

for any phase shift (p. 

Proof. For M > we compute 

M+ip M M+<p -M+ip 

-M+Cp -M M -M 

•P 

U{M + ip)- U{ip - M) dip. 







Since the integrand on the left hand side has a sign and converges pointwise, we can pass to the limit 
7\/ ^ oo by means of the Monotone Convergence Theorem, and this gives the desired result. □ 



Notice that the proof of Lemma [3T9l is very close to that of Lemma 1271] and that the assumptions on 
U can easily be weakened. For instance, ((24|) holds also for all functions U G BV as these can be written 
as the difference of two monotone functions with well defined asymptotic values. 



Lemma 13.91 implies that C is invariant under phase shifts, and hence inf £|c = inf /^|c(,. 
Corollary 3.10. For all W £ C and all phase shifts (p we have C{W{- + (p)) ~ C,{W). 
Proof. We define the monotone functions 

Ui{^) inf \W^{^), C/2M := - inf ^{{AW){^)), C/3 := - Ui, := -^AW) - U2, 

and find 

4 p 

£{Wi- + p)) - CiW) =Y.J ^^^^ + ^) - '^'^(^) = ^('J'(-l) - *(!)) = (25) 
thanks to Lemma 1321 d 



Note that the energy condition (E) is essential for Corollarv l3.10l In fact, if it is violated, then ([^5)1 
implies that C is unbounded from below and above. 

Next we derive explicit bounds for C(W) in terms of || W II2 by comparing with the functional 

C^iW) := / L{W, W)-L{Wsi,, Wsh)dip = / g^iWiip)) dip, (26) 



where the last identity holds by definition of see (|16p . and (E). In particular, (A) implies C^{W) > 0. 
Lemma 3.11. There exists a constant c < 4 such that 

\Ci^iW) - C{W)\ < c (27) 

for all W eC. 



12 



Proof. The monotonicity of W £ C implies 

W{ip-^) <W{ip) <W{(p+^), W{ip~^) < {AW){ip) <W{(p+^) 
for all (ys e R, and due to |$'(ti;)| < 1 for all w G [—1, 1] we also have 

\^{{AW){ip)) < \{AW)iip) -W{ip)\. 

In combination we obtain 

\<i>{{AW){^)) - $(iy(^)) \<W{ip + ^)- - i) 



and Lemma 13^ yields 

J \^{{AW){tp)) - ^{W{ip)) I d^ < 2. (28) 

R 

Therefore, 

\CiW) ~ C#{W)\ < J \^{iAW)icp)) ~^{Wiip)) \ dip + J |$((^W,h)M) -$(Vt^.hM)| d^<4, 

R R 

where we used (f28l) both for W and Wsh- □ 



The following lemma is key to our approach and strongly relies on the signed are condition (A) and 
the monotonicity of all profiles W £ C. 

Lemma 3.12. There exist two constants c > and c > such that 

c\\W\\l<£#{W)<c\\W\\l (29) 

for all W eCo andW Wsh - W. 

Proof. By Assumption 13.11 the hmction g$ is smooth and positive in (—1, 1) with 

5*(±1) = 0, gU±l) = 0, 5*(±1) = 1-*"(±1), 

and hence we can choose c and c such that 

c(— 1 — w)'^ < g<i>{w) < c( — 1 — w)'^ for all — 1 < w < 0, 

c{+l - wf < g^{w) <c{+l - wf for all < u) < +1. 

This implies c{W{ip))'^ < gg,{W{ip)) < c{W{(p))'^ for all LpeR, and 1^ follows thanks to 1^. □ 

Notice that the particular values of c and c depend on the choice of the pinning point (zero in Co), 
and that there is no similar result for W E C. In fact, for any W £ C we have C^{W{- + 0)) — C^[W) = 
C{W{- + (^)) - C{W) = but II W{- + Cp) II2 - II W II2-* 00 as ^ -> 00. 

Combining Lemma 13.111 with Lemma 13.121 we finally find the desired bounds for C. 

Corollary 3.13. The sub-level sets of C iti^Cq are bounded in the following sense. For each ci there 
exists C2 > such that C{W) < ci implies \\ W \\2< C2. 
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3.5 Existence of minimisers 

Now we can complete the existence proof for monotone supersonic fronts. To this end we study weakly 
convergent minimising sequences for and exploit the following observations: The limit 7 of the kinetic 
energies of the minimising sequence and the kinetic energy 7^0 of the weak limit satisfy A7 = 7 — 700 ^ 0- 
Consequently, the weak limit is a minimiser of C if and only if the limiting potential energy difference is 
less than A7. To prove this, we first observe that the potential energy converges strongly on compact 
sets so that we only have to control the potential energy in the tails of the profile. Second, we show that 
the loss of potential energy in the tails is strictly less than A7 provided the shock is supersonic. 

Theorem 3.14. L attains its minimum in C. More precisely, C is weakly lower semi- continuous on Cq, 
and each minimising sequence has a strongly convergent subsequence. 

Proof. We first prove the minimisation in C which consists of the following steps. 

1. Let {Wn)n C C be a minimising sequence, i.e. C{Wn) inf C\c. According to Corollarv l3.10l we can 
suppose Wn e Co, and Corollary 13.131 shows that the sequence {Wn)n C Co with Wn = Wsh — Wn 
is bounded in L^. Hence we can extract a (not relabelled) subsequence such that 

Wn weakly in L^ and 7„ := i |1 Wn 7 > 7oo 5 || Woe h 

for some weak limit Woo corresponding to Woa = Wsh — Woo- Moreover, we have Woo 6 Cq and 
Woo 6 Co as Co is closed under weak convergence. Setting Un ■= Wi — Woo = Wn — Woo we find 

Un-^0 weakly in L^, and ^ \\Un\\l^ J ~ loo, (30) 

which implies 

AUn and AWn AWoo strongly and pointwise in L^([-Af, M]) (31) 
for each M < 00, see Lemma [3731 With these notations we write 



C{Wn) - CiWoo) = J i((M/oo + Unf " M^^) - [^AWoo + AUn) " HAWoo)) dip 



— h, n + I2, n — h, 71, 

where the /i.„'s are given by 



12, n ■■= J WooUn - ^' (AWoo) AUn dip = J (Woo ~ A^' iAWoo))Un dp>, 

R E 

13, n-= / <^(AW^)iv){i'^Un){ip))dip, 



and $(yivFoo)(y)(') pointwise normalised potentials as in Remark 13. 51 

2. Towards an estimate for /3,„, we fix e > sufficiently small, and M = M{e) sufficiently large, such 
that 

sup \WM ~ {AWoo){^)\ < ^e. (32) 

\ip\>M 

The convergence provides AWn{±M) -> AWooi±M), hence 

liAWoo - AWn){-M)\ + liAWoo - AWn)iM)\ < \e 
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for all sufRcicntly large n, and since Wn is increasing we conclude that 



sup mh{^)-{AW„){^)\<e. (33) 

\ip\>M 

Combining ([5^ and (|33p with Remark (|3.5p we find 

< <i>Mw^)(^)((^f/„)(^)) < |(^[/„)M' (34) 
for |(^| > M and all sufficiently large n, where 

Ce := sup {$"(«;) : w £ [-1, -1 + e] U [1 - e, 
With ([M)) and \\AUn Il2<|| t^n II2 we now estimate 

0</3.n< J $(^vv^)(^)((^t/„)</?)d^+ J '^(AW^){v){{AUn)<f) dip 

\<p\<M \v\>M 

<iC2 J {AUnKv^f + ^Ce J {AUn){vfd^ 

\<p\<M \<p\>M 
<iC2 UAUn)\[-M,M]\\mCe \\Un\\l 

and exploiting (|3ip and ([50)) we obtain limsup/3_„ < Ce(7 ~ 7oo)- Since e was chosen arbitrarily, 

n — >oo 

we derive 

limsup/3,„ < Co(7 - 7oo), (35) 

n — >oo 

where Co = max{<I>"(— 1), $"(1)} satisfies < Co < 1 on account of Assumption 13. II 
3. Combining (I35p with /i,n ^ 7 — 700 and l2,n ^ we find 

lim C{Wn) > C{W^) + (1 - Co)(7 - 7oo) > C{W^), (36) 

n— *oo 

and since (W^)^ is a minimising sequence, we also have 

CiWoo) > lim /:(H/„). (37) 

n — ^00 

This shows C{Woo) = lim £(W„), so Woo is a minimiser of £ in C. 

n — >oc 

Revising the proof so far yields the additional results: (i) C is weakly lower semi-continuous on Cq as 
(|36p holds even for non-minimising sequences, (m) Concerning strongly convergent subsequences, (|36p 
and (j37p provide 7 = 700 and hence || W^„ ||2"*|| Woo II2, which in turn implies that the weak convergence 
Wn — Woo — > is strong. □ 

We have now finished the existence proof for action-minimising fronts with monotone profile. In 
particular, the second part of Theorem 11.11 follows by combining Lemma 12.21 Remark 12.31 Lemma 13.71 
and Theorem 13. 141 

We conclude with two remarks about the necessity of conditions (S) and (A). 

(i) Suppose that $ satisfies (E) with strictly concave-convex Then the asymptotic states are 
subsonic with $"(±1) > 1 and 17$ is negative in (—1, 1). Using we show that holds for 
negative constants c, c < 0, and since (|27|) still holds we conclude that C is bounded from above but 
unbounded from below Jn particular, action minimising fronts cannot exist and each minimising 
sequence diverges via || W ||2^ 00. On the other hand, similar to the proof of Theorem 13 . 141 we can 
show that C attains its maximum in C. However, we cannot conclude that the maximiser solves the 
front equation, compare Remark 13.81 and numerical simulations indicate that the maximiser is not 
a front. 



15 



(ii) More generally, suppose (A) fails because there exist —1 < w < +1 with g,s>{w) < 0, and define the 
sequence Wn G C with extending plateau at w by Wn{^) = sgmp for \ip\ > n and Wn{^) = w for 
\(p\ < n. Then and (P?)) imply C{Wn) ^ —oo, so action-minimising fronts cannot exist. Of 
course, this does not disprove the existence of local minimisers of C, but the numerical simulation 
in i j4. 21 indicate that these do not exist. 

4 Qualitative properties and numerical computation of fronts 
4.1 Qualitative properties of fronts 

In this section we discuss some aspects concerning the shape of front profiles W, which hold for any 
bounded and monotone solution W = Wsh — W to the front equation W ~ A^'{AW). Recall that such 
solutions belong to as noted in Remark 12.31 

Front profiles are actually strictly monotone for strictly convex <f> for the following reason: First we note 
that there exists at least one (fo such that W'{(po) > 0, and with W' > this implies {AW'){(p) > for 
all (p G {fQ — ^, ipo + ■^). Moreover, the front equation implies 

W = A['i>"{AW)AW'], 

and since <f>" > wc infer that W'{(p) > for all (p E {tpo — 1, cpo + l). Finally, by iterating this argument 
we find W'{ip) > for all ^ G R. 

We next study the convergence to the asymptotic states. Heuristically, the converg^cc to the constant 
states w = ±1 is dictated by the linearisation in the asymptotic states, that is W = X±A^W, with 
A+ = $"(1) and A_ = $"(—1) for right and left tails, respectively. Differentiating twice with respect to 
if we find the linear advance-delay differential equation 

W" = X±AW, 

with discrete Laplacian {AW){ip) = W{(p + 1)+W{(p — l)~2W{(p). The exponential ansatz W{(p) — e^"^"^ 
gives the characteristic equation 

= 2A±(cosh(T) - 1), 

where decay corresponds to solutions r with negative real part and monotonicity requires r G K. Since 
supersonic speed implies < A± < 1 in the normalised potential it is straightforward to see that the only 
nonnegative real roots of the characteristic equation are the double root t = and a simple root t± > 0. 

In conclusion, for monotone W the function W is expected to decay exponentially with rate t+ and 
T_ as 93 — > +00 and (p — > —00, respectively. We prove a weaker variant of this heuristic result for ip> — > 00, 
and mention that an analogous result characterizes the decay for Lp — > —00. 

Lemma 4.1. Let W be a monotone solution to pOp . Then, for all r, t with t < < t there exist 
positive constants c and c such that 

cexp {—Tip) < W{(p) < cexp (— r</?) (38) 

for all ip > 0. 

The recently established existence of centre-stable manifolds for nonlinear advance-delay differential 
equations, [Geo08| Theorem 5.1, shows that nonlinear decay is precisely described by the linear decay 
except for the centre directions. Since t+ is an algebraically simple eigenvalue there is no secular growth 
and so Lemma 14.11 in fact implies that (|38p holds for t = r = r+ and more precisely that there are 
constants a 7^ 0, > such that 

W{ip) ^ acxp{~T+ip) + C'(cxp(-(t+ + i^)ip). 
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Our point, however, is to obtain the decay rate using a direct approach with a rather simple proof. Not 
surprisingly, the same rates of decay describe the tails of small amplitude supersonic solitons in [FP99j 
Indeed, Lemma 14.11 requires monotonicity only in the tails and therefore also applies to all supersonic 
solitons with monotone tails, as these solve ([TU]) after a suitable rescaling, compare [Her09| . 

Proof. We first show the existence of both exponentially decaying lower and upper bounds for W and 
prove afterwards that r and r can be chosen arbitrary close to r+ . 

1. By Taylor-expansions of $ around w = 1 and due to {AW){(p) as ^ oo we find for each 
< (5 < 1 some (pg such that 

X+{l-5){A'^W){p>)<Wiip)<X+{l + 5){A^W){ip) for all ip>ips. (39) 

Moreover, since W is monotonically decreasing for 1^9 > we have || W \\2> Jq W{ip)'^ dip > (pW{ip)'^, 
and hence 

0<W{^) <\\W\\2^-^^^. (40) 

2. An important building block for the upper bound is the implication 

U{(p)<C^-^^^ Vv?>0 =^ {A'^U){ip + l) <C(p^'^/^- V(/9>0, (41) 

which follows from a direct computation of ^^[iy9"^/^]. For fixed S and ip > ipg wc combine the 
previous resuhs as follows: (gOl) and ([41]) imply {A'^W)((p + 1) < || h'P''^^^ and with ((391) we find 

o<w{ip + i) < x+ii + s^whip-^/"^. 

Using gl]) for U{ip) = W{ip + 1) we obtain {A^W){ip + 2) < A+(l + d)\\ W ||2^~^/^ and ^ yields 

< W{ip + 2) < Xl{l + Sf\\W\\2ip~'^^^. 
Iterating these arguments gives 

< W{(p + n)< Xlil + 5)"|| W II2 (^-^/^ V ^ > (^5 V n G N, 

and since S was arbitrary we have established the upper estimate in ()38p for some r with r > 
— In A_|- > and some sufficiently large c. 

3. Concerning a lower bound, the monotonicity of for (/? > implies {AW)['p + j) > jW{(p), and 
therefore {A'^W){(p + |) > j^W{ip). Combination with ^ gives 

w{ip + ^)>^{x+{i-5))w{ip) yip>^s, 

and by iteration we show the existence of an exponentially decaying lower bound. 

4. Using the previous step we now assume that ([38]) holds for some constants c , c, and some rates 
< T < T+ < T, and show that these rates can be improved. In relation to the above characteristic 
equation we define the function g{T) 2T~^(cosh(T) — 1) and observe that ^^[e^^*^] = p{t)c^'^'^. 
Hence, for all r > and iy9» > we have 

U{ip) <Cexp{-Tip) Vipxp^ =^ {A'^U){ip) <Cg{T)exp{-TLp) yip>ip^ + l. (42) 

Due to (|39p and using an iteration argument similar to the above one we arrive at 

W{ip) < C {X+{1 + 6) g{T)f Clip {-Tip) yip>ips + nyneN, (43) 
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and since W is decreasing this implies 

for all (5 > 0, r > 0, > and some constant C. Therefore we can improve the assumed decay rate 
r provided that A+£i(t) < 1, that means r < t+. More precisely, choosing S sufficiently small and 
adapting the constant c we derive from (|43)) that the upper estimate in ((38)) also holds for the new 
rate 

Z = T-ln(i(A+^.(r) + l)) >T. 

Since the iteration of the map t t yields a strictly increasing sequence that converges to r+, 
we conclude that r can be chosen arbitrarily close to r_|_. Finally, reversing the inequality signs on 
both sides of the implication (|42)) we can show that r can also be chosen arbitrarily close to t+ . 

□ 

4.2 Numerical computation of fronts 

In view of Theorem 13.141 it seems natural to approximate fronts by using some discrete counterpart of 
the gradient flow of C, as for instance the explicit Euler scheme ([23]) . 

A corresponding numerical scheme on a finite interval is readily derived and implemented, but from 
the rigorous point of view the account of such a scheme is limited: (i) There is no convergence proof, (ii) 
Due to the lack of uniqueness results it is not clear whether or not our existence result covers all fixed 
points of the scheme. 

Nevertheless, such schemes work very well numerically and provide moreover some intuitive under- 
standing for why energy conservation, supersonic front speed, and area condition are necessary for the 
existence of action minimising fronts. 

We start with functions <&' that have exactly one turning point. A positive example with admissible 
potential $, initial profile Wq = Wsh, and A = 1 is plotted in Figure [5] After a small number of iterations 
the profile W becomes stationary and satisfies the front equation up to high order; the same qualitative 
behaviour can be observed for other < A < 1 and other initial data Wq € C. 

In Figure[n]wc plot the result for the same scheme applied to a potential that does not satisfy the energy 
conservation. Here the profiles do not converges to a front solution, but instead form a travelling wave for 
the iteration: After some initial iterations we find a stationary profile that is shifted in each step. Recall 
that in this case the action functional C is not invariant under shifts, and thus the profiles successively 
decrease their action by converging to the asymptotic state with higher potential energy. More precisely, 
assuming T{W) = W{- + (p) and exploiting ^ we find < A/: = C{W)~C{T[W]) = (p{^{l) - «'(-l)), 
and hence sgn((^) = sgn(<I>(l) - <i>(l)) < 1. 

In order to understand the role of supersonic fronts we next consider potentials that conserve the 
energy with concave-convex Such potentials are prototypical for subsonic conservative shocks and 
have unbounded C as discussed at the end of iJH Figure [7] illustrates that the iteration scheme generates 
a plateau with increasing length and two counterpropagating travelling waves that connect the plateau to 
the asymptotic states. The height of the plateau is the unique solution to w = ^'(iD) with — 1 < w < 1, 
and as above we conclude that the decrease in C is given by A£ ~ <^i(<f>(u>) — $(—1)) 4- (^2(^(1) — '^{w)), 
where Cpi > Q and (^2 < arc the pase shifts for the travelling waves. 

Finally, we study two potentials that both have three turning points in the distance jump of supersonic 
conservative shock data. The numerical simulations in Figure [5] and Figure [5] indicate that the signed 
area condition is truly necessary for the existence of action minimising fronts. If this condition fails 
the profiles create again an extending plateau and two counter-propagating travelling waves, where the 
plateau height w satisfies ^'{w) = w and ^"{w) < 1. 
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Figure 5: Example 1 with <I> as in Assumption 13.11 The profiles converge to a fi'ont. 
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Figure 6: Example 2 with <I>(1) ^ $(—1): The profiles minimise the action through subsequent shifts and 
converge to a constant. 
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Figure 7: Example 3 with $"(±1) > 1: The profiles converge to a constant via an extending plateau. 
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Figure 8: Example 4 with area condition: The profiles converge to a front. 
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Figure 9: Example 5 without area condition: The profiles generate an extending plateau. 
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